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Determinations of \Vud\ and \Vua\ along with their implications for the unitarity test \Vud\ + 
\Vus\ + \Vub\ = 1 are discussed. The leading two loop radiative corrections to neutron /3-decay 
are evaluated and used to derive a refined relationship \Vud\ T„(i + Sqa) = 4908(4)sec. Em- 
ploying I Kid I = 0.9740(5) from superallowed nuclear decays and the measured neutron lifetime 

^N) ' r„ = 885.7(7)sec, leads to the precise prediction gA = 1.2703(8) which is compared with current 

^ I direct experimental values. Various extractions of \Vu3\ are described and updated. The long ac- 

^^ , cepted Particle Data Group value of |Kis| from fitted K^a decay rates suggests a deviation from 

04 ' CKM unitarity but it is contradicted by more recent experimental results which confirm unitarity 

CO ' with good precision. An outlook for possible future advances is given. 

^ '. 

O ■ PACS numbers: 12.15.Hh,14.20.Dh,13.20.Eb 

o 

r-| ■ I. INTRODUCTION 

O ' The study of nuclear beta decays played an important historical role in unveiling universal properties of weak charged 

D , current interactions and in helping to establish the SU{2)l x U{1)y Standard Model of electroweak unification. In 

i-G ' the limit of neglecting electroweak loop corrections, a special subset of such decays, the superallowed 0^ -^ 0+ 
Fermi transitions, depend only on the vector current which is conserved and, therefore, not renormalized by strong 

"'~j • interactions at q^ ~ 0. Hence, they are ideal for extracting \Vud\, a cornerstone of the CKM (Cabibbo-Kobayashi- 

rS , Maskawa [1113) three generation quark mixing matrix, 

Vud Vus Vub 
X^CKM = I Vcd Vcs Vcb I . (1) 

Vtd Vts Vtb 

Indeed, those decays currently provide the very precise determination Q 

\Vud\ = 0.9740(5) (0+ ^ 0+ /3-decays), (2) 

which we will discuss later in this paper. Combining the value in eq. ^ with knowledge about \Vus\ from kaon and 
Hyperon decays along with the fact that \Vub\ — 2.1(3) • lO^'^ 4] is negligibly small, allows one to confront the CKM 



unitarity relation, 



|Kd|' + |Ks|' + |K6|' = l, (3) 



at a high precision level. 

That prediction has been tested to about ±0.15%, an impressive accomplishment. At that level, it has confirmed 
the presence of very large, ^ 4%, Standard Model loop corrections in the extraction of \Vud\ from the data [a,|a]. 
However, for many years, a small deviation, ^ 2a, from exact unitarity has been persistently observed. That issue has 
been further clouded by various distinct \Vus\ determinations which seem to be inconsistent with one another. These 



problems are addressed later in some detail, when new K^z results fully consistent with unitarity are also discussed. 
Of course, if a real deviation from unitarity expectations is seen, it signals the presence of as yet unaccounted for new 
physics beyond the Standard Model, an exciting prospect. Alternatively, if unitarity is respected, constraints on new 
physics can be implied. However, \Vud\ and |T/us| must be thoroughly scrutinized both theoretically and experimentally 
before conclusions are drawn. 

Neutron beta decay, n — > pe^V^, has not until recently been prominent in efforts to determine \Vud\ and test 
unitarity. It depends on both vector and axial-vector charged current interactions. The latter are renormalized by 
strong interactions at q^ ~ 0. The size of that effect is parameterized by gA = Ga/Gv, a fundamental quantity in its 
own right. Indeed, the value of qa, which has grown over time from about 1.23 to 1.27, is important for predicting 
the expected solar neutrino flux Q , light element abundances from primordial nucleosynthesis [g , the spin content of 
nucleons [^ [i3| and for testing the Goldberger-Treiman relation [IJl . A byproduct of the analysis in this paper will 
be to provide a very precise determination of qa that can be compared with more direct neutron decay asymmetry 
measurements of that important parameter or to refine the above mentioned applications. 

As new more intense neutron facilities turn on, experimental measurements of both the neutron lifetime r„ and gA 
from the electron asymmetry in polarized neutron beta decay are expected to become much more precise T3|. Indeed, 
combining determinations of those two quantities, can yield \Vud\ with an anticipated uncertainty competitive with 
the error in eq. |(5J), i.e. dominated by theory. In preparation for those improvements, we present in this paper a 
relationship among \Vud\, Tn and gA which includes one and some dominant two loop quantum corrections. It can be 
used to determine gA from t„ and the \Vud\ input from superallowcd nuclear /^-decays or as an independent measure 
of \Vud\ using T„ and neutron decay asymmetry determinations of gA- 

Our plan is as follows: in section ^ we update the radiative corrections to neutron decay by incorporating the 
©(a^) effects due to leading logs, (some) small next-to-leading logs and Coulombic effects. The last of these has 
been considered previously, but with the wrong sign. We take this opportunity to correct that longstanding error. In 
section UTTl we review and update (slightly) the extraction of |Kid| from superallowed /3-decays. Using that value of 
\Vud\ along with the neutron lifetime r„ we derive in section HVI a very precise prediction for gA and compare it with 
direct asymmetry measurements of that parameter. Then, in section^ we review and update various determinations 
of \Vus\ and point out inconsistencies among them. The main problem stems from old K^s decay rates obtained from 
fitted PDG studies and is suggestive of errors in some longstanding (accepted) kaon decay properties. Indeed, recent 
Kes results from Brookhaven and Fermilab experiments confirm significant errors in the old charged and neutral kaon 
decay branching ratios and lead to values of \Vus\ fully consistent with unitarity. Implications of a unitarity violation 
or confirmation in eq. Q are briefly discussed and an outlook for future advances is given in section IVTI 

II. RADIATIVE CORRECTIONS TO NEUTRON DECAY 

Our analysis of the radiative corrections to neutron beta decay builds on the results of earlier studies, particularly 
the classic work by Wilkinson [1^. They included 0{a) radiative corrections as well as effects due to the final state 
electromagnetic ep interaction embodied in the Fermi function. A number of other small corrections from proton 
recoil, finite nucleon size etc. have also been examined |14L llSl lla | . 

In the Standard Model, one renormalizes the beta decay amplitude using p^ lla. Il9l l20l l2lL 12^ 

G^ = 1.16637(1) • 10"''^ GeV~^ (4) 

the Fermi constant as obtained from the muon total decay rate (inverse lifetime) . In that way, ultraviolet divergences as 
well as radiative corrections common to both decay amplitudes are absorbed into G^. The remaining loop differences 
and bremsstrahlung effects can then be factorized into an overall 1-l-RC correction to the neutron lifetime, r„ 

1 _ Gl\Vy,d\ 5 . 2 



T„ 27r^ 

where / is a phase space factor, 



-m 



{l + 3gl)il + RC)f, (5) 



/ = 1.6887, (6) 



which includes a relatively large Fermi function contribution |l3l | (~ 5.6%) as well as smaller nucleon mass, size and 
recoil corrections. It has been somewhat updated in eq. © to incorporate slight nucleon mass shifts. 

We note that the electroweak radiative corrections, denoted by H-RC, have been factorized in the same way for 
both the vector and axial-vector contributions. (Interference and induced coupling corrections are negligibly small in 
the case of the lifetime [iJiliS-) That factorization effectively defines gA via the relative normalization of the axial- 
vector current as measured by the lifetime (it incorporates QED as well as strong interaction effects in its definition) . 



Employing such a definition for gA means tliat tfiere will be some 0{a) corrections to the g^^^ measured in neutron 
decay asymmetries that must be applied before contact with eq. (O can be made at a level of high precision. (Those 
corrections will be discussed in section Hvl 1 In this connection, it is worthwhile to note that the normalization of 
gA in eq. Q is consistent with the so-called 1/fc method, in which the radiative corrections to various observables 
in /3-decay, such as the lifetime, the electron spectrum, and the longitudinal electron polarization are expressed in 
terms of effective couplings Gv and G a [ ^ [241 l2fiL l2fl . The same approach has been employed to calculate the 
corrections to the electron asymmetry zr, /^ l29l l30l Isij . The factorization of the short distance contributions of 
the radiative corrections, implicit in eq. |(SJ, conforms also with an asymptotic theorem concerning their behavior in 
arbitrary semileptonic decays mediated by W^ |33| ■ 

To order a, the radiative corrections in eq. JSJ are given by [sL la. 1331 13^ 

RC = ^ |5(i?„0 + 41n^+ In ^ + 2C + A,1, (7) 

zvr 1^ TOp niA J 

where g{Em) represents long distance loop corrections and bremsstrahlung effects averaged over the /3-decay spectrum 
[23i] {Em = 1.292579 MeV, the end-point electron energy in neutron decay), 

^9{Em) = 0.015056. (8) 

The other parameters in that expression have the values 

mz = 91.1875 GeV, 



TOp = 


= 0.9383 GeV, 


niA : 


- 1.2 GeV, 


C c 


- 0.891, 


A, c 


- -0.34. 



(9) 

Here, following 0j|2a,l23) we have approximately identified niA with the mass of the ai(1260) axial vector meson. 
The value of C is an update of the calculation discussed in Ref. ^^ (using gA = 1-27 for self consistency). This leads 
to a first order result 

RC (©(a)) = 0.03770, (10) 

which is rather sizeable. (Together with the Fermi function, such (primarily) QED corrections increase the neutron 
decay rate by 9.37%!) 

Given the magnitude of the order a corrections in eq. IjlOf) and our desire for high precision, it becomes imperative 
to include the leading 0(a^) contributions and estimate the theoretical uncertainty in the radiative corrections. 
Regarding the latter issue, the overall uncertainty is usually obtained by allowing uia in eq. (U)) to vary by a factor 
of 2 up or down. That reflects the fact that the last 3 terms in eq. result from axial-vector current loop effects 
and their calculation is not perfectly matched in going from long to short-distance contributions. The scale, niA, 
uncertainty reflects the matching error in a rough, but numerically realistic, way. With that methodology, the 
theoretical uncertainty is estimated to be 

RC (uncertainty) = ±0.0008, (11) 

an error common to all neutron and nuclear /5-decay studies. Reducing that dominant theory matching error would 
be very useful, but it is extremely challenging and beyond the scope of this paper. 

Our focus in this section is to include the so called leading log corrections of the form a" In" ^^ and a" In" ^^ , 

n — 2,3, . . . in eq. ^ along with some of the other potentially most important Ola^ln^^j and 0[a^) effects. 

Regarding the last of those, there is a relatively important Coulomb correction not included in the product of the 
Fermi function and 1-l-RC. It is called ^<5 in the literature and approximated by (see Appendix) 

—S:::i-a^ln^- + ... = -OmOAi. (12) 

27r TOg 

Although considered previously, for some reason it was given the wrong sign. Correcting the value -)-0.0004 used in 
the past to eq. (|12() corresponds to a —0.00083 shift which is quite significant at the level of our analysis. 



In the Appendix, we give formulas that sum the leading log contributions in eq. (Uj) via the method in ref. [35 
They lead to the replacements 

1 + — In— ^ S{mp,mz) = 1.02248, 



TT m. 



■p 
1 + T^ln— ^ -> L(2^„,mp) = 1.02094. (13) 

where the large |^ In ^^ contribution is hidden in the ■^'g{Em) function of eq. Q f23'|. 

Our final new input is to estimate next-to-leading log (NLL) corrections of the form a^ In ^-^^ and a^ In ^^p coming 
from fermion vacuum polarization insertions in loops with photon propagators. Because they all enter with the same 
sign and there are quite a few leptons and quarks that contribute, one might expect those fermion loops to dominate 
the NLL contribution. However, as illustrated in the Appendix, they turn out to be quite small in the MS formalism 
we employ. We estimate 

NLL = -0.0001. (14) 

Other C(q!^) contributions are not expected to be significant, but a complete calculation to C(q;^) would be very 
difhcult and beyond the scope of this paper. Also, such a refined calculation is not obviously warranted until the one 
loop matching uncertainty in eq. Hll|l is significantly reduced. 

Our next step is to organize the radiative corrections into a factorized form that does not induce spurious contribu- 
tions when multiplied out. To accomplish that end, we take an effective field theory approach, dividing contributions 
into very long distance (the Fermi function), long distance, intermediate distance, and short distance factors. Of 
course, we must make certain that the matching is done correctly, for example by introducing terms such as -^S (see 
eq. C2J) when appropriate. Overall, we employ the following factorization beyond the Fermi function 



1 + RC = 



L{2E^,mp) + -C+—S 

TT ZTT . 

S{mp,mz) + ^^(ln'-^+Ag]+NLL\, (15) 

2t: \ TTiA J 

where a{mp) ~ 1/134. Employing the values we have given above for the quantities in eq. p5|l . one finds 

1 + RC = 1.0390(8), (16) 

where the uncertainty in eq. Hll() has been employed. Comparing eqs. (|16() and H10|) . one sees that the C'(a^) and 
summation effects have increased the RC by 0.0013, not a very large shift. Nevertheless, they must be included in 
precision studies. In summary, our updated analysis includes the following two-loop, C(a^) contributions: 1) Leading 
Logs, 2) Next- To-Leading Logs from fermion vacuum polarization effects and 3) Coulombic matching corrections due 
to the factorization of the Fermi function. Other neglected two-loop corrections are not enhanced by relatively large 
factors and thus assumed to be negligible. 

Employing eq. I|l()|) in eq. ©, we derive the relationship 

|Kd|'(l + 3.gi)r„ = 4908±4sec. (17) 

That master formula can be used to extract \Vud\ via 



/ 4908(4) sec 



1/2 



as T„ and gA become more precise or to obtain gA from t„ and |V^d|. Currently, one finds from the experimental 
averages 

<™ = 885.7(7) sec (19) 

5^1™ = 1.2720(18), (20) 

the CKM parameter 

|Kd|=0.9729(4)(ll)(4), (21) 



where the errors correspond to r„, gA and RC respectively. In section lllll we compare that value with the more precise 
\Vij,d\ obtained from superallowed beta decays. Of course, the full utility of eq. (|18|) will be much better realized as r„ 
and gA measurements improve. 

The master formula in eq. 117|l can also be used to extract a very precise value of gA via 



9A 



1636(1) sec 1 



(22) 



Note that the theory uncertainty (±1 sec) in eq. (|22|) is due to the niA scale uncertainty in the radiative corrections 
and it cancels with the same error in determinations of \Vud\ from nuclear decays. The formula in eq. H22f) will be 
utilized in section Hvl 

III. SUPERALLOWED /3-DECAYS AND |Kd| 

Superallowed 0+ -^ 0+ Fermi transitions have been the focus of many studies, both experimentally and theoretically. 
Here, we start with the recent results of Towner and Hardy l36l |. They have thoroughly scrutinized the nine very 
well measured superallowed /3-decays, using RC last updated in ref. [3^ [S^, HI], and taking great care to correct 
for various nuclear Coulombic and structure dependent effects. In the end, they arrive at J-t values that are nucleus 
independent and can be used to extract |l/ud|, modulo uncertainties in the radiative corrections. The Z independence 
of their results is an important consistency check, since the daughter nuclei have Z values ranging from 5 to 26 with 
correspondingly different Coulomb and structure corrections. 

In Table m we give the (slightly) updated values of \Vud\ obtained from the nine best measured superallowed /3- 
decays, incorporating the isospin symmetry-breaking and structure dependent corrections Sc and Sns used by Towner 
and Hardy, in conjunction with our new factorization scheme for the radiative corrections in eq. (|15|l (changing Em 
and S in that expression as appropriate for each nucleus). Specifically, we include the corrections 6ns in the second 
factor of eq. (|15|l and append an additional factor (1 — 5c). 

TABLE I: The values of \Vud\ obtained from superallowed beta-decays. The error given does not include nuclear and theory 
uncertainties common to the analysis. It is taken from ref. 36] and used to obtain the weighted average. 



Nucleus 


IKdl 


loc 


0.97388(76) 


14 Q 


0.97445(41) 


26^1 


0.97416(35) 


3*C1 


0.97431(40) 


38^ 


0.97424(43) 


*2Sc 


0.97351(38) 


46 y 


0.97372(43) 


^^Mn 


0.97396(44) 


5*Co 


0.97409(43) 



The values of \Vud\ derived from those distinct measurements are very consistent with one another and range from 
about 0.9735 to 0.9745. The weighted average is centered at \Vud\ = 0.974047; so, we round down to 

|Kd| = 0.9740(l)(3)(4), (23) 

where the uncertainties are experimental, nuclear theory and RC (see cq. (|11|) '). We have checked that the combined 
small changes in the RC arising from our new factorization in eq. (|15() . improvements in the higher order leading 
logs, central value of tua, small NLL (—0.0001) correction etc. tend to cancel. For that reason, the result in eq. H23|) 
is essentially the same as the one obtained earlier by Towner and Hardy (cf. eq. (j^J) using the radiative corrections 
given in Ref. [33 ■ We note that the sign of the ^(5 corrections for the superallowed decays (which are all e+ emitters) 
is correct in the literature and numerically the change in sign of that correction for neutron decay was our biggest 
modification of previous results. So, \Vud\ in eq. (|^ remains the current best value. Indeed, comparison with \Vud\ 
extracted from neutron decay t„ and gA values in eq. (|21l) shows that they share the same RC uncertainty, but the 
error on gA must be improved by about a factor of 4 before the neutron decay becomes competitive. 



One way of reducing the RC uncertainty in \Vud\, perhaps by as much as a factor of 2 would be to use the 
TT"*" —^ TT^e'^Uf. decay rate for which the loop induced axial- vector contributions are better controlled and nuclear 
theory uncertainties are circumvented j^. However, the small branching ratio ~ 10~® makes that method statistically 
challenging. Nevertheless, an ongoing PSI experiment finds 



|Kd|= 0.9749(26) ■ 



BR(7^^ 



-^e(7)) 



1.2352-10-4 



1/2 



(24) 



where its dependence on the tt"^ 



Veij) branching ratio (used for normalization) is exhibited. We assume the SM 



theory value of 1.2352 • 10"'* is correct and hence obtain \Vud\ = 0.9749(26), in excellent agreement with nuclear and 
neutron results but with a larger error. Alternatively, others have chosen to use the PDG recommended branching 
ratio of 1.230(4) • 10~* which leads to \Vud\ = 0.9728(30) which is also in agreement within errors even though the 
central value may appear low |40| . 

Note, if we average the \Vud\ determinations above, we still find 

iKdl = 0.9740(5), (Average) (25) 

because the superallowed /3-decays dominate and they have an average central value slightly larger than 0.9740. 

IV. qa: theory VS. EXPERIMENT 



Employing the neutron lifetime measurement in eq. 119|l and the value of \Vud\ in eq. (|23|) . we find via eq. (|22|l the 
Standard Model prediction 



gA = 1.2703(6)(5), 



(26) 



where the errors are in t„ and \Vud\ (nuclear uncertainty). That precise prediction is larger than the PDG recommended 
value of 1.2670(30) |41,], but smaller than the single best asymmetry value 42] of 



gA = 1.2739(19), 



(27) 



(after small QED corrections are applied). 

For now, the qa in eq. (|26l) should be considered the standard, to be used in solar neutrino flux calculations etc. 
However, to take proper advantage of its precision, one should employ the neutron decay definition we have assumed 
and correct the weak interaction process under consideration for its own electroweak radiative corrections. As an 
example, consider the relationship between the qa defined via the neutron lifetime whose value is given in eq. H26() 
and one defined by the lowest order polarized neutron decay asymmetry 



A = 



2g(i-g ) 

1 + 3.9 



?,2 



'^ asy 

9 = 9a ■ 



(28) 



If we wish to replace g by gA in that expression, then there will be additional (energy dependent) radiative corrections 
to the asymmetry. Those corrections were computed long ago, originally by Shann and later confirmed by many 
others. They are quite small, leading to about a 0.1% shift in the asymmetry |23|- That effect is corrected for in 
the most recent [425 experimental measurement of the asymmetry, although it is well below current experimental 
uncertainties. It would be useful to compute the QED corrections to other processes where the very precise value of 
gA in eq. (|26|l might prove useful, for example primordial nucleosynthesis or the solar neutrino flux. 



Wusl DETERMINATIONS AND CKM UNITARITY 



The 2002 PDG recommended value for |Ks| Eu, 

|Ks|pDG = 0.2196(26) PDG2002, 



(29) 



has remained fixed (modulo small variations in its uncertainties) for many years. It is based on fitted branching ratios 



from rather old data for K^s decays, K^ 



andK+ 



Tj^+ . The resulting K^s decay rates are proportional to | Vu. 
analogous to superallowed nuclear beta decays (or tt^ 



n"e+iy, 



> TT e+i/e, combined with the kaon lifetimes, tkj^ and 
and can be used for its extraction. In fact, they are 
= ) in that only the weak vector current contributes at 



the tree level. Since that current is conserved in the SU(3) flavor limit, strong interaction corrections are of second 
order in SU(3) breaking. Those effects, characterized by the departure of the form factor /+(0) from 1 along with 
isospin breaking effects were considered in the classic study by Leutwyler and Roos 43] that forms the basis for the 
extracted value of |T4s| in eq. H29|l . 

In addition to SU(3) breaking effects, there is a fairly significant first order rrid — triu correction due to n'^-ri mixing 
(~ 4%) that must be separately applied to the charged K^^ decay rate and an extra Coulombic 7r^e+ final state QED 
interaction for K^^. The fact that even with those different isospin violating corrections, both the neutral and charged 
Ke3 decay rates gave consistent values for \Vus\ (eq. (|29|l is their average) has often been used to argue for the validity 
of eq. H29|l as compared with for example Hyperon beta decays which have tended to give somewhat larger values for 
\Vus\ but are not as theoretically clean. 

We note that combining eq. if^ with the value of \Vud\ in eq. (P5|l and using \Vub\ — 2.1 x 10^^ leads to 

|Kdlo+^o+ + lK.|pDG + \Vnbf = 0.9969 ± 0.0015. (30) 

That roughly 2 sigma deviation from unitarity has been a persistent problem for many years. It has been at times 
interpreted as a hint of new physics or as an indication that something is wrong with the data and/or theory 
calculations used to extract \Vud\ and \Vus\- Since perfect unitarity would require a rather large shift (+ 3.2%) 
in \Vus\ but a relatively smaller shift (+ 0.16%) in \Vud\, the latter has been generally thought to be the root of the 
problem. As a result, considerable experimental and theoretical scrutiny has been applied to |V^d|. Nevertheless, as 
emphasized in the first part of this paper, its value has remained rather stable. 

To illustrate the above approach and some of its underlying uncertainties, we describe the general K^j decay rate 
formula Hal, 



1927r3 



r{K ^ nei^i-f)) - ^^^^EW (1 + SI,) C' |Ks|' /^(O)/^, (31) 



where C^ = 1 for Kl or Ks decays (to both ir^e^) and C^ = 1/2 for K^. Sew — 1.022 is the universal short-distance 
radiative correction in eq. (|15|l [32l | while 6]^ are model dependent long-distance QED corrections recently estimated 
to be (for radiative inclusive studies) jij; HE ^M 

Sy = +1.3 ±0.3%, 

S<j^+ = -0.1 ±0.7%. (32) 

The form factor /+(0) incorporates SU(3) breaking. Leutwyler and Roos found 

/+(0) = !+"''' = 0.961 ± 0.008, (33) 

a value recently confirmed by a lattice calculation l47l an d to some extent by new papers based on Chiral Perturbation 
Theory (ChPT) pa. lia. |49|. In particular, in Ref. |48| it is shown that the only unknown constants in the ©(p^) 
contributions in ChPT can be determined by accurate measurements of the slope and curvature of the scalar Ki^ 
form factor. In Ref. ^0] dispersion relations are employed to calculate theoretically these two observables, leading to 
/+(0) = 0.974 ± 0.0057 ± 0.0028 ± 0.009, which is consistent with Eq. ((23, although the central value and estimated 
errors are somewhat larger. In this paper we employ the classical result given in Eq. (|33|l . but we also emphasize the 
importance of refining the lattice and ChPT calculations of /+(0). 

In the case of charged kaons, md — ruu mass splittings give rise to tt — 77 mixing such that 

jK+^'' ~ 1.022/:[^"'^" = 0.9821 ± 0.008 ± 0.002. (34) 

Finally, the phase space factor is determined for a linear form factor to be (for a slope factor of 0.028) 

4,0 = 0.1550, 

/^+ = 0.1594, (35) 

while for a quadratic form factor suggested by KTeV data [51| , 

/|,o = 0.1535. (36) 

The difference, ~ 1%, is somewhat accounted for by assigning an extra ±0.7% form factor decay rate uncertainty. 
However, we note that the central value of \Vus\ will depend on which parametrization is used. The value of \Vus\ in 
eq. (|29ll would become 0.2207 if a quadratic parametrization were employed. 
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Recently, a number of developments have cast doubt on the reliability of eq. (|29|) . First, a new measurement of the 
K^^ branching ratio by the E865 Collaboration at Brookhaven National Laboratory finds it to be about 5.3% larger 
than the fitted PDG value used to obtain \Vus\- Their finding conflicts with the earlier K^s decay rates, even after all 
isospin violating differences are taken into account and on its own leads to |52j 

|Ks|=0.2236(23)//:f^^'^"(0), E865 (37) 

where f^ "^ (0) = 1+SU(3) breaking effects, including m^ — m^ corrections (but unlike E865, we have not absorbed 

QED corrections into its definition). For the value f!^ ^ (0) — 0.9842(84) effectively used by the E865 collaboration 
after removing the QED correction, one finds jV^^I — 0.2272(30) and correspondingly almost perfect unitarity. 



|Kdlo+^o+ + IKsIe865 + IK6|' = 1.0003 ±0.0017. (38) 

However, that is not the end of the story. In an even more recent development, the KTeV Collaboration (E832) |53 
at Fermilab has reported a thorough analysis of all primary Kl decay modes, including Kg^. They find significant 
disagreement with the PDG fit values for several of the main Kl branching ratios, including K^^, exactly the types 
of shifts required to bring the Kl system into accord with unitarity and the K^ results of E865. 

The KTeV Collaboration finds (using the PDG Kl lifetime) 

r{KL -^ neiy) = 0.520(4) x lO^^-* MeV. (39) 

That radiatively inclusive K^s partial decay rate is about 5% larger than the 2002 PDG fit value. After accounting 
for measured form factor effects, a new calculation of QED radiative corrections ,44 ;45:. i46j and SU(3) breaking (via 
Leutwyler and Roos |43|), they find 

|K.|= 0.2253(23) KTeVife3, (40) 

where (because of the very high statistics) the error is essentially dominated by SU(3) breaking, form factor shape 
and Kl lifetime uncertainties. Unlike K'^^, the Kl extraction is not directly sensitive to the up-down mass difference. 
For Kf^3, they obtained a similar result \Vus\ = 0.2250(23), which provides strong confirmation. Taken together with 
the value of \Vud\ in eq. lf2!^ . one finds from eq. (|^ 

\Vudf + IK.I' + iKftl' = 0.9994(14) KTeV K,3 , (41) 

in very good agreement with unitarity, a remarkable turn of events. So, E865 and KTeV are in apparent agreement 
regarding \Vus\ and unitarity. Of course, they are dependent on the K'^ and Kl lifetimes which, if the history of the 
kaon branching ratios is any indication, could change upon closer scrutiny. 

It is useful to average the K^s and Kf^^ results from E865 and KTeV in order to get a single determination of 
\Vus\ and unitarity constraint that can be used to limit new physics appendages to the standard model. To carry 
out such an average requires a consistent treatment of features common to the analyses of both measurements. In 
that category we put SU(3) breaking effects and the form factor shape (linear vs. quadratic). We adopt the KTeV 
quadratic form factor along with its ±0.7% uncertainty and assume the Leutwyler-Roos estimate of SU(3) breaking 
and jTid — nriu effects. Together they shift the E865 value of |Kjs | up by about 0.7%. So, we wind up with the following 
two quantities to be averaged: 

IK^I = 0.2288(26)(20) Shifted E865 A'+g, 

iKsl = 0.2252(13)(20) KTeV /f°3 and ^Og average, (42) 

where we employ the KTeV Kg.^ and iC^a averages under the assumption of muon-electron universality. The common 
error from the form factor and SU(3) breaking, (±20), has been factorized. The remaining uncertainty in the K^ 
case comes from adding in quadrature the experimental uncertainties along with errors due to QED, the K^ lifetime, 
md—mu effects, and normalization uncertainties associated with correlations among K^ branching ratios. We estimate 
the last of these to be ±0.5%. For the Kl, the first error comes from experimental uncertainties, the Kl lifetime 
error (which is appreciable), and QED. 

Carrying out a weighted average, using the first set of errors to do the weighting, we find 

iKsl =0.2259(12)(20) E865-KTeV Average, (43) 

or, combining errors in quadrature, 

iKsl =0.2259(23) E865-KTeV Average. (44) 



Together with the \Vud\ value in eq. (|2l, that gives 

iKdl' + IK.I' + iKbl' = 0.9997(10)(10), (45) 

where the first and second errors arise from \Vud\ and \Vus\, respectively. Combining the errors in quadrature, we 
obtain 

iKdl' + iKsl' + iKbl' = 0.9997(14), (46) 

which confirms unitarity superbly. Agreement with the standard model expectation can be used to constrain various 
new physics effects. For example, if the muon had some additional, exotic decay mode such as ^ ^ e+ the wrong 
neutrinos, it would affect unitarity through the G^ used in our normalization pj]. The good agreement in eq. H46(l 
limits the branching ratios for those types of hypothetical decays to be < 0.2%, which is similar to the best direct 
constraints, but more general. 

At this point we note that the central values of \Vus\ in eq. H42|l differ by about 1.2a. It could be a simple effect 
due to experimental systematics (statistical errors are very small) or might indicate a problem in the Kaon lifetimes, 
rud — TTiu, QED long distance radiative corrections or K'^ correlated branching ratios. Improved measurements in the 
charged and neutral kaon properties may prove useful in clarifying the difference. However, the total uncertainty in 
\Vus\ of about ±1% is dominated by SU(3) breaking and the form factor shape and magnitude. It will be difficult to 
reduce those errors much further. Unitarity prevails, but unless some new procedure for calculating SU(3) breaking 
effects with much higher precision is found, it seems unlikely that the error in eq. (|44|l can be reduced much further, 
i.e. a ±1% uncertainty in \Vus\ is near the end of the road for K^^. Similar remarks apply to \Vud\ where theory 
uncertainty dominates. Fortunately, they seem to be ending with a triumph for unitarity and a strong confirmation 
of the standard model. 

Other new analyses of i^es by the KLOE Collaboration at Frascati and NA48 at CERN are being completed and 
should report results shortly. It will be interesting to see if they confirm E865 and KTeV or reopen the K'es problem. 

In another relatively recent development, Cabibbo, Swallow and Winston (CSW) |53,|53 have revisited the extrac- 
tion of \Vus\ from Hyperon beta decays. That procedure is sometimes criticized as unreliable because the decay rates 
are renormalized by first order SU(3) breaking effects in the axial-vector contributions. However, rather than just 
employing total decay rates, CSW also included decay asymmetry measurements that effectively measure the first 
order SU(3) breaking effects, analogous to the use oi gj^^ in neutron beta decay to determine |Kid|- They found 

iKs I =0.2250(27) Hyperon Decays, (47) 

where the error quoted is purely experimental and SU(3) breaking has been neglected. Nevertheless, if taken at face 
value it gives 

|Kd|o+^o+ + IK^Inyperon + IK^l' = 0.9993(16), (48) 

i.e. good agreement with unitarity. Because SU(3) breaking effects and various other theory uncertainties have not 
been considered, we do not include eq. (I47|l in our averaging. 

Finally, a recent lattice calculation j^] of the pseudoscalar decay constants 

Jk/U = 1.210(4),tat(13)syst, (49) 

can be combined [Sg with the experimental quantity 



to yield 



1.3336(44), (50) 



Wu. 



\Vud 

where the uncertainty is lattice dominated. If unitarity is assumed, that ratio implies 



0.2278(26), (51) 



iKsl = 0.2221(24), 

iKdl = 0.9750(5), (52) 



10 

while if \Vud\ = 0.9740(5) is employed 

|Ks| =0.2219(25) Lattice. (53) 

It corresponds to a modest 1.4 sigma deviation from unitarity. The beauty of the lattice approach is that it is still in 
its infancy. With new larger computers and better treatment of chiral symmetry and isospin breaking, the uncertainty 
in eq. H53|l may be reduced by a factor of 4 or more, making the lattice approach to /k/ /tt ultimately the best way to 
determine \Vus\- Should that happen, it will be somewhat ironic that the axial current determination via K12 and 7:12 
decay constants turns out to be theoretically more pristine that the vector current approach using Kes- If something 
similar were possible for \Vud\, the unitarity test (constraint) in eq. H46|l might be significantly improved. 

Recently, hadronic r decays have been employed to determine \Vus\- The value found in |59l |. using LEP data, is 

iKs I = 0.2208(34) T decays. (54) 

This determination will likely improve with new data coming from BaBar and Belle. 

The \Vus\ central values discussed above vary from about 0.22-0.23 depending on the data used and SU(3) breaking 
corrections applied. That range is to be compared with the value suggested by 0+ — > 0+ beta decays and perfect 
unitarity 

IK^I =0.2265(22) Unitarity + 0+ ^ 0+ Nuclear. (55) 

A comparison of the different determinations of |y„s| is illustrated in Fig. ^ 



PDG 2002 ' 

Shifted PDG 2002 ^ 
Tau lepton ' — 
Lattice 



Hyperon > • > 

KTeV I • 

Shifted E865 ' • ' 

Unitarity and Nuclear beta decays > • > 

0.21 0.215 0.22 0.225 0.23 

FIG. 1: Determinations of \Vua\ from various sources. The Hyperon value does not include theory errors. Shifted values 
correspond mainly to the change of linear to quadratic form factor parametrization. All K13 results assume Eqs. 133I34II . 



VI. CONCLUSION AND OUTLOOK 

We have presented an update of the electroweak radiative corrections to neutron /?-decay. It currently provides 
the most precise determination of gA, a quantity that finds many applications in nuclear, particle, and astro-physics. 
As QA and t„ experimental measurements improve, our results can be used to obtain a \Vud\ that is competitive 
with nuclear /3-decay determinations (which yield |T4d| = 0.9740(5)), but without the nuclear structure dependent 
uncertainties. In the end, we will still be limited by a ±0.0004 theory error that comes from uncertainties in the 
axial- vector induced loop corrections. Reducing the latter error further will require a new theoretical approach or 
the measurement of F(7r+ — > Tr^e^j^e) which has smaller uncertainties. The 10~^ branching ratio makes the latter 
strategy very difhcult statistically. 

More problematic than \Vud\ m testing CKM unitarity has been inconsistencies in the various \Vus\ determinations. 
Of particular concern has been the low value ^ 0.22 obtained from PDG fits to old Ke3 data which for a long time 
has suggested a small but persistent departure from unitarity. However, recently measured Kes branching ratios 
for both the charged and neutral kaon exhibit large deviations (about 5%) from the PDG fit values, increasing 
\Vus\ to a level consistent with unitarity. In fact, the average value they provide |T4s|=0. 2259(23) together with 
|Ktd|=0. 9740(5) from superallowed beta decays concur with unitarity expectations up to ±0.0014 (see eq. (|46ll '). That 
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good agreement confirms the predicted large +4% [3 radiative correction of the Standard Model at the 28a level! 
A triumph for quantum loop effects. However, before one can be sure that the unitarity problem has been resolved, 
the recent results should be confirmed by anticipated new measurements of the K^s branching ratios from the KLOE 
experiment at Frascati and the NA48 experiment at CERN. In addition, given the recent changes in kaon branching 
fractions, it would be nice to have new confirmation measurements of both the K^ and K^ lifetimes. 

The newest approach to \Vus\ appears to be very promising |^. It combines a lattice calculation of Jk/ f-K with 

the experimental measurement of pV ^^ ^ ^ , J . The latter ratio is already very well measured, but in view of 

the changes in K^ decay rates, the numerator should be experimentally revisited. Also, the clectroweak radiative 
corrections that largely cancel in the ratio should be reexamined. 

The lattice calculation of Jk/ f-n has some very nice features that suggest its already small uncertainty may be 
reduced much further. They include a cancellation of the statistical uncertainties which are highly correlated and the 
scale uncertainty. The challenge is to refine extrapolations to the chiral and continuum limits with full dynamical 
fermion simulations. That should become possible with the advent of very large special purpose computing facilities. 

The current confirmation of unitarity embodied in eq. H46(l gives us further confidence in the Standard Model and 
can be used to constrain "New Physics" appendages to it. Roughly speaking, it rules out additional tree or loop 
level contributions to the \Vud\ + |Kts| combination at about the ±0.2% level. Effects of about that magnitude 
could in principle come about from supersymmetry loops ,6^ l6lL l6^ , an extra Z' boson [63] , hea vy q uark or lepton 
mixing, exotic muon decays (since we normalize in terms of G^), compositeness, extra dimensions [gJ etc. However, 
there are no really compelling reasons to expect such a large deviation. For that reason, confirmation of unitarity via 
experimental or theoretical changes has been anticipated. Its verification is, nevertheless, important. Also, finding 
new ways to improve the unitarity constraint, e.g. lattice calculations, is difficult, but should be strongly encouraged. 
Such challenges push our computational and experimental skills to new limits, stimulate ingenuity, and perhaps at 
some point new physics will be uncovered. 

Note added: After the completion of this paper, a new low temperature measurement of the neutron lifetime was 
reported, r„ = 878.5(7) (3) sec [63. It differs from the world average in Eq. ifT^ by about 6.7cr. In conjunction 
with \Vud\ = 0.9740(5), it would lead on its own to a larger gA value, namely gA = 1.2766(6)(5); instead, if used in 
combination with the best asymmetry value gA = 1.2739(19), it would imply \Vud\ = 0.9757(4)(11)(4). Although this 
value is not in sharp disagreement with Eq. (j^J (the difference is about 1.3(t), clarification of the neutron lifetime 
differences by new improved experiments is clearly an important goal for the future. 

Also, the KLOE collaboration has presented some preliminary new results on K,.^ and K^-^ decay rates and the Kl 
lifetime |7(ll | . They confirm the branching ratio increases observed by KTeV (within errors) and when finalized should 
improve the Kl lifetime world average. 

Finally, a recent analysis JZII by the NA48 Collaboration at CERN finds a neutral Kez decay rate consistent with 
the KTeV result in eq. (|39|l . However, they effectively employ (after extrac ting electromagnetic contributions) a larger 
value of /+(0) = 0.974 suggested by recent chiral perturbation results |4^ Hq, 113, IH, H^, 113 and a linear form factor 
parametrization of phase space. As a result, they obtain |l/us| = 0.2187(28) which suggests an overall 2.2a deviation 
from unitarity. Such a significant change in the interpretation reinforces the importance of refining lattice and chiral 
perturbation theory calculations of /+ (0) as well as the need to better experimentally determine its q^ dependence. 

APPENDIX 

In this appendix, we describe the input that went into: 

1. the leading log summations L{2Ejn,mp) and S{mp,mz)] 

2. the partial next to leading log calculation due to photonic vacuum polarization insertions and 

3. the residual -^5 correction that results from a proper matching of the Fermi function and long distance 0{a) 
corrections. 

1. Leading Log Summations 

We follow the approach of ref. ,33j where a renormalization group summation for the leading short-distance logs was 
given and a value for S{mp,mz) derived. Here we extend that method down to the intermediate region 2£',„ — irip. 
That simply requires a change in the anomalous dimension from 2a/ n to 3a/27r and an evaluation of the MS (modified 
minimal subtraction) coupling a(/Lt) at low scales. 
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The leading log summation in an MS approach is simply given by 

L(2i?™,m,) = f^^Y^' (ap)Y^'' fap^Y^'' (a(jn^Y^'' (^X^'\ (56) 



a{2Em)J \a{mu) J \a{md) J \a{m^)J \a{ms) ^ 

/I ^\3/4 / / N\9/16 / / \\9/19 / / \\9/20 / / n x 36/17 

^K,-z) = f^) f^) f^) f^) f^) , (57) 

Va(mp)/ Va(TOc)/ V"("^r)/ \a[mb)) \a{mw) J 

where the a(/i) are values of the MS (modified minimal subtraction) QED coupling at a scale ^. This coupling is 
given to leading log order by 

a-i(M) = a-\m,)-^y2Q}e{^l^mf)h-i^ + ^e{^i-mw)ln-^, (58) 

dTT ^ — ' ■' rnr Ztt mw 

/ ^ 

a-\me) = 137.036 +-!- = 137.089. (59) 

In that expression, the sum is over all quarks and lepton flavors / (with a color factor of 3 for quarks). We thereby 
find 

a-\2E„^ = 2.585 MeV) = 136.745, 

a-i(m„ =62 MeV) = 136.0708, 

a"i(md = 83MeV) = 135.9263, 

a-^{m^) = 135.7896, 

a-i(TO^ = 215MeV) = 135.2368, 

a-^{mp) = 133.9861, 

a-\mc = 1.35 GeV) = 133.67728, 

a-^irrir) = 133.3662, 

a"i(TO6 = 4.5GeV) = 132.1174, 

a-i(TOH' = 80.4 GeV) = 128.0389, 

a-i(TO2 = 91.1875 GeV) = 128.001. (60) 

Relatively low effective quark masses 1^3 have been used in those results in order to incorporate QCD contributions 
at low energies. In that way a~^{mz) = 128.0 is obtained while a more detailed higher order QED and QCD analysis 
including e^e^ — > hadrons via a dispersion relation gives a^^{mz) = 127.934. 
Employing the above values in eq. H56() and eq. H57() . we find 

L{2E,n,mp) = 1.02094, 
S{mp,mz) = 1.02248. (61) 

Those results are not very sensitive to the quark masses employed. For example, changing the to„ value by a factor 
of 2 leads to a shift in the RG by less than 3 x 10~^ which is well below the overall ±8 x 10"^ uncertainty assumed 
in the RC. 

A useful relation, valid for the nine transitions in Table QJis 

L(2i?^,mp) = L(me,mp)(^^^|^j , (62) 

which leads to 

L(2Ern, rrip) = 1.026725 1 1 - ^^^ In ^^ . (63) 

V Stt rrie J 

2. Next-to-leading logarithmic corrections due to photonic vacuum polarization. 

Next-to-leading logarithmic (NLL) corrections to weak decays have been studied in detail for QCD (see, for example, 
|6a.l67l |l. Here we adapt those results to QED. We consider a subset of contributions, containing fermion loop insertions 



13 






FIG. 2: Vacuum polarization diagrams contributing logarithmic and NLL corrections to the beta decay rate. QED loop 
corrections to external charged particle legs are included but not illustrated. 

in the photon propagator, as depicted in Fig. 13 That particular subset is somewhat enhanced by the large number 
of contributing fermions. Therefore, we expect it to dominate the NLL contributions. 

The effect of short distance corrections can be parameterized as follows. We keep only leading-logarithmic terms, 
and those NLL ones that arise from fermion loops. Then, the decay width F can be expressed in terms of the 
lowest-order width F° times a correction factor, 



F 



pU 



1 



■7 



(°)La- 



= ln^ 



La^ (7(1 
47r ' 



2poB 



1 



-L^a^^^o) 



(/3o + 7 



(0) 



(64) 



and we will be interested in the value of this correction at fi : 



A similar analysis can be extended down to 2Er, 



The leading order anomalous dimension of the four- fermion operator m7''(1 — 75)^(8)67^(1 — 75);^ is given by diagrams 
similar to those in Fig. |21 but without fermion loops. It can be expressed in terms of fermion charges. 
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(0) 



[Qu - Qdf + QeiQe + 8Qu - 2Qd) 



(65) 



The running of the coupling constant is described by a sum over all contributing fermions with charges Qf and number 
of color varieties Uf (equal 3 for quarks and 1 for others). 



Po 



■3"' 



/ 



(66) 



At the NLL order, we also need the two-loop anomalous dimension 7'^' and the matching coefficient B. The former 
is determined by calculating 1/e poles of the diagrams in Fig. [3 The latter is given by the finite part of analogous 
diagrams without fermion loop insertions. Their individual values depend on the scheme of the calculation (for 
example, treatment of 75), and only their combination occurring in eq. (|64|l is scheme-independent. 

We list here the values of both quantities obtained in the naive dimensional regularization (NDR), with an anti- 
commuting 75, and in the 't Hooft-Veltman scheme (HV), 



^Sl 



7ndr 



B 



HV 



23 
6 ' 



B 



NDR 



19 
6 ■ 



In both schemes we get 



2(3qB + 7 



(1) 
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With these results, the logarithmic corrections of eq. (|64|l become 

Mw 



F 



F" 1 



2Q(/i) ^^^ Mw 



-) In; 



M 



^ + 2 In 



M, 



w 



M 



4n 



(67) 



(68) 



(69) 



All but the last terms in this correction factor are included and summed up to all orders in the factor 5* in eq. H57|l . 
The last term is the fermionic NLL correction. Because it is very small, we estimate it as though all fermions u, d, s, c, b 
and e, /z, r contributed over the whole range from Mw to mp. In this approximation, we have 



20 



(70) 
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and the numerical value of the NLL fermionic correction is 

4n 



4n /aN 



\ Mw 
m 



-0.00007. 



(71) 



Small additional NLL contributions from the Em to rap region shift that value to about —0.0001. Other two-loop 



NLL and non-logarithmically enhanced corrections are individually of order 
we assume that they can be neglected. 



■In 



2E„ 



0.00003 or smaller. Hence, 



3. Evaluation of the 0(Za^j corrections 

In this Appendix we summarize the evaluation of the O(Za^) corrections in the case of neutron /3-decay (Z = 1). 
These are defined as the residual corrections of this order not contained in the product of the Fermi function and the 
0{a) corrections. For a point nucleus, analytic results for the correction to the charged lepton spectrum has been 
obtained in the extreme relativistic approximation |6S| |. The sign is opposite for electron and positron emitters (as 
can be readily understood by a glance at the figures in Ref. JGg) and, for the neutron decay, we have 



AP 



In 



43 5 /2E 
18 " 3 ^ VtoT 



(72) 



(see Eq. (5) of Ref. JGgl). It has also been verified that, to good approximation, the leading contributions to Eq. I|72p 
extrapolate smoothly to their non-relativistic limit, which has also been obtained analytically. Next, we evaluate 
numerically the spectral average of ln(2Ty) {W = E/m) using the phase space factor \/W'^ — 1 W (WVax — W) , 
leading to 



(ln(2W)) = 1.1390. 



Combining Eqs. (|72|l and 173|l . we obtain 

(AP> 



-8.006a^ 



-4.26 X 10' 



(73) 



(74) 



In order to take into account the finite proton size, we employ Eqs. (7-9) of Ref. [33, with the sign again reversed, 
since we are dealing with an electron emitter. These formulae correspond to a uniformly charged sphere of radius 

R ~ vfaj where a = V6/A is the rms charge radius of the proton. Using a = 0.90 fm [l^, we have A/nip — 0.572. 

Insertion of this value in Eqs. (7-9) of Ref. J33|, leads to a finite proton size effect —5.4 x 10~®. Combining this result 
with Eq. (|74|l . we obtain our final answer for this class of corrections in the case of neutron decay: 



a 
— ( 

2tt 



-4.3 X 10""*. 



(75) 



We have verified that terms of Ol (A/m^) ], not contained in Eqs. (7-9) of Ref. |33, as well as estimates of the 
corrections OyZ^a^j, give negligible contributions in the case of neutron decay. 
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